Effective Gauge Field Theory of the t-J Model in the 
Charge-Spin Separated State and its Transport Properties 
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We study the slave-boson t-J model of cuprates with high superconducting transition temperatures, 
and derive its low-energy effective field theory for the charge-spin separated state in a self-consistent 
manner. The phase degrees of freedom of the mean field for hoppings of holons and spinous can be 
regarded as a U(l) gauge field, Ai. The charge-spin separation occurs below certain temperature, 
Tcss, as a deconfinement phenomenon of the dynamics of Ai. Below certain temperature rsG(< 
less), the spin-gap phase develops as the Higgs phase of the gauge-field dynamics, and Ai acquires a 
mass mA. The effective field theory near TsG takes the form of Ginzburg-Landau theory of a complex 
scalar field A coupled with Ai, where A represents d-wave pairings of spinous. Three dimensionality 
of the system is crucial to realize a phase transition at TsG- By using this field theory, we calculate 
the dc resistivity p. At T > Tsg, p is proportional to T. At T < Tsg, it deviates downward from the 
T-linear behavior as p oc r{l - c(rsG - TY}. When the system is near (but not) two dimensional, 
due to the compactness of the phase of the field A, the exponent d deviates from its mean-field 
value 1/2 and becomes a nonuniversal quantity which depends on temperature and doping. This 
significantly improves the comparison with the experimental data. 
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I. INTRODUCTION 

In many physical systems, one may identify their 
microscopic Hamiltonians, but calculations of physical 
quantities starting from these Hamiltonians may not be 
straightforward. One promising way to overcome this dif- 
ficulty is to derive an effective theory that is appropriate 
for the energy scale in question and use it for calculations. 

The discovery of cuprates with high superconducting 
transition temperatures Tc has stimulated condensed- 
matter physicists for one and a half decades. Here the t-J 
model serves as a canonical model of cuprates, and is ex- 
pected to have potentiality to explain their various inter- 
esting properties observed in experiments. The t-J model 
itself can be regarded as a low-energy effective model of 
the Hubbard model or the d-p model at the energy scales 
of the hopping amplitude of electrons i(~ 0.3eV), and 
the antiferromagnetic (AF) spin couplings J(~ O.leV), 
which are much lower than the energy scales of the strong 
onsite Coulomb repulsion [/(~ 3eV). The t-J model al- 
lows us to understand some basic physical properties of 
strongly-correlated electrons. For example, the origin 
of the superconductivity is attributed to the formation 
and condensation of hole pairs at nearest-neighbor (NN, 
sites in the background of short-range AF spin order 
Furthermore, the mean-field theory (MFT) in the slave- 
boson (SB) representation predicts an interesting phase 
diagram in the doping ((5)-temperature (T) plane such as 
the spin-gap state in which the NN AF spin pairings de- 



velop. However, it is still not easy to calculate physical 
quantities such as dc resistivity analytically from the t-J 
model itself. Thus it is welcome to derive an effective 
theory of the t-J model at lower energies. It is preferable 
that this effective theory takes a form of local field theory 
so that many established techniques are applicable. 

Many quasi-two-dimensional cuprates exhibit certain 
anomalous metallic behavior above Tc in various quan- 
tities such as-| dc resistivity. Hall coefficient, magnetic 
susceptibility.0 For example, the dc resistivity p at fixed 
5 shows the T-linear dependence on T above certain 
temperature. These anomalies call for a new theoretical 
explanation, probably in a framework beyond the conven- 
tional Fermi-liquid theory. Andersono pointed out that 
the charge-spin separation (CSS) phenomenon may ex- 
plain them. 

The CSS is a phenomenon in which charge and spin de- 
grees of freedom of strongly-correlated electrons behave 
independently. In the SB representation, an electron is 
viewed as a composite of two constituents, a holon and a 
spinon, which bear charge and spin degrees of freedom of 
electron respectively. The CSS is "naturally" described 
by the SB (or slave- fermion) MFT, in which the holons 
and spinous are treated as quasifree particles having no 
correlations among them. When one incorporates fluctu- 
ations of phase degrees of freedom of MFs, these phases 
behave as gauge fields coupled to holons and spinous, and 
the system possesses a lo&al U(l) gauge symmetry. 

In the previous papersQ (hereafter we call them Pa- 
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pers I), we have studied the possibihty of the-CSS in 
the t-J model using gauge-theoretical methodsE By in- 
troducing auxiliary fields, the system can be viewed as 
a lattice gauge model. The corresponding gauge dynam- 
ics has two general possibilities in its realization; (i) the 
confining phase in which holons and spinous are confined 
to electrons and the gauge fields fluctuate strongly, or 
(ii) the deconfining phase in which holons and spinous 
appear as almost independent quasi-free particles and 
the gauge-field fluctuations are weak and can be treated 
by the usual perturbation theory. We interpreted the 
CSS as the second possibility above; CSS is a deconfine- 
ment phenomenon of gauge theory. Furthermore, we ob- 
tained the result that the system exhibits a confinement- 
deconfinement phase transition (CDPT) at certain tran- 
sition temperature Tcss i below which the deconfinement 
phase takes place, i.e., the CSS takes place below Tcss- 
The CDPT is of second order when the system has three- 
dimensional couplings whatever small they are, while it 
is infinite order of Kosterlitz-Thouless type when the sys- 
tem is purely two dimensional. 

Very recently, study of the CSS was revived, and 
the confinement-deconfinement problem of gauge theories 
of nonrelativistic fermions is addressed. Unfortunately, 
most of these studies do not quote the previous works in 
which important results on that problem were obtained. 
Moreover, some of the recent arguments are in apparent 
contradictions to the previous results, but these authors 
do not discuss the origins of these contradictions. For 
example, in his recent papeiO Nayak discussed the above 
problem and concluded that the slave particles (holons 
and spinous) are always confined by U(l) gauge field. 
This result is obviously in sharp contrast to our result in 
Papers I which predicts a deconfinement phase at low T . 

How does this contradiction come out? In Nayak's 
paper, dynamics of the gauge field is not discussed at all 
and it is simply concluded that the infinitely strong gauge 
coupling in the bare Lagrangian necessarily leads to con- 
finement. On the other hand, in Papers I, we studied the 
gauge-field dynamics by using nonperurbative methods. 
Even if the gauge coupling is infinite in the original lattice 
model, the gauge field can acquire nontrivial dynamics at 
low energies such as a deconfinement phenomenon due to 
the couplings to matter fields. In short, due to ample pair 
creations of holons and/or spinous, fluxes of electric- like 
field connecting external charges which would cause the 
linear-rising confining potential are truncated into short 
segments, giving rise to a deconfining potential. In lattice 
gauge theory such as multi-flavor QCD,eI this fact is now 
well-established and verified by numerical studies. We 
showed that a similar deconfinement phenomenon occurs 
in the gauge theory of the slave-particle t-J|jnodel. 

Abcwe is the main point of our commenta to Nayak's 
paperlj .Another point we argued is that the Elitzur's 
theoremtJ does not exclude the possibility of deconfine- 
ment phase. There are some misunderstandings that the 
Elitzur's theorem prohibits the existence of the deconfine- 
ment phase, since the theorem states that the averages 



of gauge-noninvariant quantities should vanish. However, 
MFT, for example, can describe the deconfinement phase 
in accord with the Elizur's theorem by pjeraging over the 
gauge-rotated pspies of a MF solution.!!^ 

In his reply|13 Nayak admitted that these two points, 
are certainly correct. There he also posed an argumentEHI 
that there are ambiguities in assigning EM charges to 
slave particles, which makes the possibility of CSS quite 
doubtful. This is a misunderstanding. Our explanation 
is as follows; This ambiguity is related to the choice of 
reference state from which the EM charges are measured. 
If measured from the vacuum, the EM charge of a holon is 
zero and the charge of a spinon is the same as the charge 
of an electron, e(< 0). Similarly, if they are measured 
from the half-filled state, the holon charge is — e and the 
spinon charge is zero. Since this is an important point 
to understand the slave particle approactw-sve present our 
explanation in Sect. 6 of the present papert^l after deriving 
the relevant formulae in Sect.SA. 

In the present paper, we focus on the CSS state in the 
SB t-J model and derive the low-energy effective field 
theory. By using the lattice model suggested by this field 
theory, we confirm that the nature of the transition into 
the spin-gap state is not a crossover but a genuine phase 
transition for the three-dimensional system. We apply 
this effective field theory to |Ga,lculate the dc psistivity. 
Previously, Lee and NagaosallJ and loffc et al.ll3 showed 
that p (X T for fermions and bosons interacting with a 
massless gauge field. Their system has some relation to 
the uniform RVB MFT of the SB f-J model. Recent ex- 
periments on YBaCuO by Ito et alS3 and others reported 
that p deviates downward from the T-linear behavior be- 
low certain temperature Tsq(> Tc). This TgQ coincides 
with the temperature determined by NMR and neutron 
experiments at which a spin gap starts to develop.EZi So 
it is quite interesting to derive an effective theory in the 
spin-gap state of the SB t-J model, and calculate the 
spin-gap effect on p as a function of 6 and T. The ef- 
fective theory used in Ref.Efl is inadequate for this pur- 
pose, since it assumes no spin gap. In Sect. 3, we shall 
see that the spin gap generates a mass of gauge field via 
the Anderson-Higgs mechanism. Thus the effect of scat- 
terings of holons and spinous by gauge fields is reduced. 
This gives rise to a downward deviation of p from the 
T-linear behavior. Our explicit pssult of p in Sect.SC is 
consistent with the experimentsJlj 

The strucutre of the present paper is as follows; In 
Sect. 2, we start with the MFT of the SB t-J model, and 
review the results of Papers I for the U(l) gauge dy- 
namics of the t-J model in a compact but selfcontained 
manner. Various analytical and numerical results which 
are relevant for the present system are consulted in or- 
der to obtain the correct phase diagram. When T is 
higher than the CDPT temerature Tcssli^) which de- 
pends on the holon density 6, T > Tcssli^): the gauge 
dynamics is realized in the confinemnet phase and the 
holons and spinous are confined to electrons. When 
T < Tcss{S), the gauge dynamics is realized in the 
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deconfinement phase and the CSS takes place as men- 
tioned. The MFT predicts another phase transition at 
Tsg(< Tqss), below which the spin gap develops. This 
transition is sometimes claimed to be a crossover rather 
than a genuine phase transition.lla We shall discuss this 
problem in Sect. 4 in detail, and show that it is a genuine 
phase transition when the system has weak but finite 
three dimensionality. 

In Sect. 3, we shall derive the low-energy and low- 
temperature effective field theory in the CSS state of 
the SB t-J model. The result of Sect. 2 assures that this 
derivation is self consistent, since the gauge dynamics 
is deconfining at T < Tcss. The effective theory is a 
Ginzburg-Landau (GL) theory L{X,Ai) coupled with a 
gauge field Ai , where the complex scalar A represents the 
d-wave spinon pairing and Ai is the phase of the hopping 
amplitujde of holons and spinous. Halperin, Lubensky, 
and Mallj considered a similar system, where Ai corre- 
sponds to the electromagnetic (EM) field and A to the 
Cooper-pair field. They calculated the effect of Ai on 
A to conclude that it converts the second-order phase 
transition to a first-order one. Then it was pointed out 
that the radial fiuctuations of A may keep-the transition 
second-order in some parameter region.cJ This is sup- 
ported Iw, numerical simulations. More recently, Ubbens 
and LeeEll calculated the one-loop effect of Ai in the SB 
MFT of the t-J model, and concluded again that the 
pairing transition at TgQ becomes first order. (How- 
ever, their Tgc appears below the superconducting tran- 
sition temperature Tc, so they concluded that the spin- 
gap phase is completely destroyed by gauge-field fluctu- 
ations.) In the present study, we take the compactness 
of the phase of A into account. Even in the CSS state, 
the compactness generates nontrivial intera ction vertices 
that are missing in the previous treatments. EjO We find 
that the periodic interaction above stabilizes the system 
at TsG, which is higher than Tc at low S. This allows us 
to use this effective theory to study the spin-gap state at 
Tc < T < TsG- We find that the critical exponents (such 
as d that characterizes the gauge field mass mA through 
niA oc (TsG — T)'^) become T and ^-dependent, that is, 
they are not universal anymore. 

In Sect. 4, we study the phase transition into the spin- 
gap state for the 3D system in detail by using the lattice 
Abelian Higgs model which is a natural extension of the 
effective field theory of Sect. 3 to a lattice model. Be- 
low certain on-set temperature Tsg (which is shown to 
be lower than Tcss)5 the spin-gap phase develops as a 
Higgs phase of the gauge-field dynamics, in which the 
gauge field acquires a finite mass m^. jE.ecent numerical 
studies on 3D U(l) gauge Higgs modelE3 show a genuine 
phase transition and existence of the Higgs phase which 
supports our discussion on the spin-gap phase. In the 
region Tsg < T < Tqss, the gauge dynamics is realized 
in the so-called Coulomb phase in which tua = 0. At the 
end of Sect. 4, we also comment on the recent dual vortex 
theory of bospsiized fermions in pure two dimensions by 
Balents et al.c3. 



In Sect. 5, we calculate the dc resistivity p as an exam- 
ple of calculations of physical quantities in the spin-gap 
state by using the effective field theory obtained in Sect. 3. 
First, we obtain the expression of p in the random-phase 
approximation (RPA), which depends crucially on the T- 
dependence of niA- Then we substitute niA calculated in 
Sect. 3 to this p. 

The results on these problems of Sect. 3-5 have been 
ported in part briefiy in our previous two short letters.cil 
In Sect. 3-5 of the present paper, we account for the de- 
tails of these results in a self-contained manner. 

Sect. 6 is devoted to discussion. We present a couple 
of comments related to the problem of observing holons 
and spinons in the CSS state and ambiguity of charge as- 
signement for holon and spinon. These comments,and re- 
marks should clarify certain misunderstandingsll3 on the 
confinement-deconfinement problem in the slave-particle 
studies of strongly-correlated electron systems. 

In Appendix, we consider the U(l) Higgs model in the 
3D continuum that is tightly related to the lattice model 
we studied in Sect. 4. We study its phase structure in 
terms of dual vaciables, and critisize the recent result of 
Nagaosa and LeeEa on this model. 

II. PHASE STRUCTURE OF THE GAUGE 
DYNAMICS 

In Scct.2A, we review the MFT of the SB t-J model 
and show the phase diagram in the S-T plane. In Sect.2B, 
we study the effect of fluctuations of MFs in a general 
viewpoint of gauge theory. A CDPT takes place at Tcss, 
below which the CSS occurs. In Sect.2C, we briefly ex- 
plain the spin gap state as a "dual" Higgs mechanism of 
lattice gauge theory with multiple gauge fields. 

A. MFT of the SB t-J Model 

Let us start with the t-J modelil Hamiltonian, 

x,i 

where Cxa is the electron annihilation operator at the 
site X of a 2D or 3D lattice with the z-component of 
spin a(=t, i)- We use i{— 1, 2(, 3)) both as the dipjction 
index and as the unit vector in the i-th direction.EII The 
SB representation of C^a is written as 
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where bx is the bosonic annihilation operator of holon, 
and fxa- is the fermionic annihilation operator of spinon. 
They satisiy the local constraint 



y ] fxa-fx 



1, 



(2.3) 



so that only the three states for each x, |vac), C|||vac), 
Cj^Jvac), (Ca;cr|vac) = 0) are allowed as physical states 

and the double-occupancy state C'j^iCj! Jvac) is excluded 
due to the large on-site Coulomb energy U. These three 
states are expressed in the SB representation as bl\0), 
/itlO)' = = 0)' respectively. The 

Hamiltonian (^]^) is rewritten for the physical states in 
terms of b^, fxa- as 



H ^ -tJ2(bi^Jl^+^.bx+i^■C. 



ft ft _ 



xlJ x+il 



(2.4) 



This H preserves the total number of holons and spinous 
at each x, the L.H.S. of (2.3), and so maps a physical state 
to another physical state. In path integral formalism, the 
partition function Z{/3) = Trexp(-/3i/) [/3 = (/cbT)"!] 
is given by 



Z = J [db][df][da] exp(-S'), 



dr 



in 



bi 



dbx 
dr 



t ^ 



i?LC ^i^axibibx + ^fLfxa - 1), 



H + Hlc\ , 
(2.5) 



where r (0 < r < /3) is the imaginary time and [db] = 
Y\x dbi(T)db x{T) , etc. Here bxir) (we omitted the 
argument r in (|2.5| )) is regarded as a complex variable 
at X and r, and fx<j{T) is an anticommuting Grassmann 
variable. iJLC respects the constraint (2^)|-wa the inte- 
gration over the Lagrange multiplier ax{T)E3 

To set up the MFT, we introduce two complex auxil- 
iary fields Xxi and X.^ on the link {x,x + i) to decouple 
both t and J term^3 as 



i?MF — 



x.i 

a 

- \{>^^^{fU^U^i - fyU^^) +H.C.}]. (2.6) 

In the MFT, the local constraint ( ^.3| ) is relaxed to the 
following global constraints for averages; 



(2.7) 



where 5 [€ (0, 1)] is the doping parameter. This modifi- 
cation is vahdated when the system is in the CSS state 
a posteriori because the CSS implies that the local con- 
straint becomes irrelevant for holons and spinous. This 
is the subject of Sect.2B, and some supplementary dis- 
cussion is also given in Sect. 6. Then, Z{(3) is written 
as 



Z= I [db][df][dx][d\] exp(-5'), 
5'= / dr|>7fo,t^+E// 



Hmf + 

U 1 / 

^/^ = - E (/^^^i^- + E /-/-O ' (2.8) 

x a 

where we replaced -ffLC by g^, /ib .f in which are the 
chemical potentials to enforce ( |2.7| ). By differentiating 
the integrand in ( |2.6| ) w.r.t. Xxi and Xxi, we obtain the 
following relations (Schwinger-Dyson equations) among 
averages; 

(xL) ~ (E/ <yfxa + 'Tj^x+i^x) , 



3J 



i^xi) ~ 9 ifxlfx+il fxl-fx+il)' 



(2.9) 



which show that Xxi describes hoppings of holons and 
spinous, while Xxi describes the resonating-valence-boud 
(RVB) (NN singlet spin-pair) amplitude. 

The MFT can be set up first by parametrizing the 
auxiliary fields as 

Xxiir) = XiUxiir), Uxi{t) = exp(iyl„(T)), 

A„(t) = X^Vx^{T), ^^(r) = exp(iB„(T)). (2.10) 

by ignoring the site- and r-dependence of their ampli- 
tudes, and then ignoring the phase fluctuations by set- 
ting Axi = 0, Bxi = 0. The effects of these phases are 
considered in Section 2B. Numerical studies were per- 
formed for this MFT in 2D with various patterns of the 
MF's Xi and A,;. One of the most interesting case is the 
so called unform RVB, in which Xi are uniform, Xi = Xi 
and Xi are the d-wave configuration, Ai = — A2 = A. In 
Sect.3C, we shall study the general pattern of Xi in terms 
of the GL theory and show the stability of this d-wave 
configuration. The MF phase diagram in 5-T plane can 
be calculated in a straightforward manner. In Fig.|l| we 
show the result of Ref.El. There appears several critical 
temperatures. Below T^, x develops, x is estimated at 
small ^'s as 



at TsG < T where A:^* = 00 From (|j|) and (Q, one 
expects holons and spinous may hop independently for 
X 7^ 0. Thus, in the level of MFT, the CSS state is re- 
alized below T^. However, as we shall see in Sect.2B, 
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the effect of phase fluctuations of Xxi reduces this CSS 
onset temperature down to Tcss which is much lower 
than T^. Below rsG7 ^ develops, so the system enters 
into the spin-gap state. There is another transition tem- 
perature Tbc below which bose condensation of holons, 
{bx) 7^ 0, takes place. To obtain a nonvanishing Tbc > 0, 
a weak-three-dimensional coupling = ax is included 
in Ref.E3. In the SB MFT, the superconducting phase 
is deatribed by the simultaneous condensations of A and 
(6a;) ,e3 because the order parameter of superconductivity, 
the pairing amplitude of NN hole states, is expressed as 
follows; 

{CxlCx+i^ — Cx]Cx+ii ) 

^ {■fl]fl+il ^ flifl+i^) ^ {bxbx+i) 

= Hhx?. (2.12) 

Thus Tc{5) is the lower temperature of T\{5) and T-Qci^)- 



B. Beyond MFT: Phase Fluctuations and U(l) 
Gauge Theory 

In any MFT, it is indispensable to examine the sta- 
bility of its solution. Especially, the SB t-J model ( |2.5| ) 
possesses a time-dependent U(l) local gauge symmetry. 



hx{T) ^ e'^-(^)6,(T), 

fx.{T) e^'-<^^^fx4r) 
axir) axir) + drOxir). 



(2.13) 



The Lagrange multiplier ax (t) is regarded as a time com- 
ponent of the gauge field. 

In the partition function (|2.8| ) of the decoupled MFT, 
we fixed the gauge to the temporal gauge by setting 
(^x{t) =0. A careful and—precise treatment of this 
gauge fixing is given in RefE3, which assures us to red- 
erives the results of Papers I. After this gauge fixing, 
there still remains in the decoupled Hamiltonian (2.6) 
a time-independent residual gauge symmetry (with r- 
independent 9x); 



bxir) ^ e^«^6,(T), 

fxair) ~* e^'Vxa(T), 

X«(r) -> e-^'-Xx^{T)e^'''^ 

\x^{T)^e''^\x^{T)e''^^'. 



(2.14) 



The last two relations, which are naturally understood 
when one recalls ( p.9[ ), show that the phase degrees of 
freedom of the auxiliary fields Xxi and Xxi can be re- 
garded as two kinds of gauge fields; 



Axiir) Axiir) - 9x + Ox+i, 
Bxi{T) Bxiir) + 9x + Ox+i, 



(2.15) 



MFT, Axi = Bxi = 0, reduces to the dynamics of the 
resultant gauge theory of Axi and Bxi- 

Let us first consider the region T > Tsg, where Ai = 
and one needs to consider only Axi{T) or Uxiir). The 
region T < Tsg, where both Axi(T) and Bxiij) appear, 
shall be considered in Sect.2C. The dynamics of Uxiir) 
is governed by an effective action that is obtained by in- 
tegrating out the holon and spinon variables. This idea 
is similar to that for the 0{N) nonlinear cr-model, in 
which an effective action or a potential of the auxiliary 
(Lagrange multiplier) field is calculated by integrating 
out the original scalar fields Uxa satisfying the local con- 
straint ^xanxa = 1- The successive large N analy- 
sis fot 3D system reveals the exitence of a phase in which 
the local constraint is irrelevant. Tbie phase corresponds 
to the CSS phase of the t-J model.ES What is important 
in integrating over bx and fxa is to intact the dynam- 
ics of Uxi- For example, if one reduces the lattice sys- 
tem into a simple continuum field theory at this stage, 
the action may contain only quadratic terms of Axi such 
as the kinetic term, \drb{x)\'^ and \{di — iAi)b{x)\'^ , and 
the Maxwell term, {diAj — djAiY . This action fits well 
to a perturbative analysis in which Ai is assumed as a 
small quantity, but it necessarily loses the possibility to 
describe the confinement phase in which Ai is large, al- 
though the integration over b[x) and fa{x) can be done 
exactly. As we proposed in Papers I, we can perfom the 
integration over bx and fxa by an approximate but a 
nonperturbative method, i.e., the hopping expansion in 
powers of Uxi, which gives rise to an action that is ca- 
pable to describe both confinement and deconfinement 
phases. This action Sjj generally contains two different 
kinds of terms; the electric one Sjje that controls the r- 
dependence of Uxiir), and the magnetic one Sum that 
controls the spatial- variation of Uxiir). Typical forms of 
them are given as 

Su = Sue + Sum, 
1 /"^ 

^UE = ^rPr{^ dTdrUli{T)drUxi{T), 

9e\^ ) „ ^0 



Sum — — y 
9 m 



where Jlpi ^^2:1 denotes the product of four UxiS along 
an oriented plaquette; UxiUx+i,jUlj^^ Jjly gE{T) and 
9m{T) are the T-dependent effective coupling constants. 
In Papers I, their T-dependences are calculated explicitly 



g|(T) oc T^, gl;{T) oc const, 



(2.17) 



by using the hopping expansion w.r.t. Uxi^ 

In the conventional lattice gauge theory,E3 the action 
in the temporal gauge Axq = is given by 



both of which are associated with the common U(l) 
gauge symmetry. Thus, the problem of stability around 



cLGT _ cLGT , cLGT 
D — '-'_B "T ^M 1 
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^ JO 

XI 

pi ^° pi 

= -^E / dT\{cosF.,,,{T), (2.18) 
^ „i "'0 „, 



where F^y(T)(= A^t + A^+i j-A^+j i-A^j) is the mag- 
netic field penetrating the plaquette. In the Hamiltonian 
formaUsm,LJ the corrsponding Hamiltonian H^'^"^ is de- 
rived as 



LGT _ E/LGT I ^LGT 



E 



M J 



(2.19) 



x,i<j 



Here, the electric field Exi{= drAxi) is the variable 
canonically conjugate to Axi, so the following uncertain- 
ity principle holds; 



AE„ A Ax i > 1- 



(2.20) 



Since Axi{(z [0, 27r]) is an angle variable, E^i takes an 
integer in its diagonal representation. When one applies 
Uxi to a state, a segment of electric flux is created on 
(x, x + i); Uxi is a creation operator of electric flux. The 
time-independent residual gauge symmetry restricts the 
physical states so as to satisfy the Gauss's law. 



'^^ViExi = "^^yEx+i,! — Exi) = 0, 



(2.21) 



where is the lattice difference operator; Vifx = fx^ 

fx- 



The system ( |2.18| - |2.2lD is known to exhibit a CDET at 
finite T. This was sho wn b y Polyakov and SusskincS by 
mapping the system (2.19) to a XY spin model. Monte 
Carlo simulations confirm the CDPT. Below we present 
a general explanation why the CDPT takes place. For 
this pur pose, let us first ignore the magnetic term Hj^^"^ 
m ( p9| ), the reason of which shall be explained later. 
The partition function is then written as 



Zlgt 



OQ 

n E n'^E.v.i...,o e-p(-/3e'E^'«)- 

xi Exi — — oo ^ 



(2.22) 



For large couphng and/or at low T such as » 1, 
configurations with small \Exi\ are favored. To inves- 
tigate the gauge dynamics, let us introduce two test 
charges with Q — ±e at xi,X2- Then the Gauss's law 
is modified to ViExi — Qx = sSx,xi — £Sx.x2j ^^id the 



configuration with the lowest energy is the one with elec- 
tric fiux Exi = ±1 formed along the straight line con- 
necting xi and X2, and Exi — elsewhere. See Fig]|a. 
The potential energy of this configuration is linear rising, 
V{xi, X2) = e'^\xi — X2\, so the two charges are confined. 
I n th is confinement state, AExi = 0, so the relation 
( ^.20| ) implies AAxi = 00, that is, Axi fluctutes violently. 
On the other hand, for small coupling and/or at high 
T such as e^(3 <C 1, Exi fluctuates violently, AExi ^ 1, 
which implies AAxi ^ 1- Thus, the fluctuations of Axi is 
small, and can be treated by usual perturbation theory. 
The potential energy can be evaluated via exchange of a 
massless gauge boson, leading to the Coulomb potential. 



V{xi,X2) 



(?\x\- 



'X2\ 



in 3D space. This is the decon- 



finement state. The electric fluxes are attached isotropi- 
cally to each external charge. See Fig.^D. These consid- 
erations lead us to a certain finite temperature Tcdpt at 
which the CDPT takes place. For fixed e^, the string ten- 
sion vanishes and charged particles are deconfined above 
Tcdpt- The quark-gluon plasma in QCD at high T is 
such a state. Below Tcdpt, particles with nonvanishing 
charges that couple to Axi are confined to form charge- 
neutral objects. Mesons and baryons in QCD are such 
objects. 

Let us estimate the effect of the magnetic term H^p'^ . 
Let us consider a state in the confinement phase, Exi of 
which takes definite values for all links. When H^p"^ 
is applied to this state, each segment of electric flux 
at (x, x + i) is deformed to the one lying on a detour 
(x, x + j),{x+j,x + j + i),{x + j + i,x + i). Thus H^p^ 
lets Exi to fluctuate and increases AExi, and so acts 
to favor the deconflnement phase. Thus Hj^'^ may de- 
creases Tcdpt certainly, but Tcdpt does not vanish, as 
the strong coupling expansion w.r.t. ^ 1 at 11,= 
assures us the existence of the confinement phase. Ej 

Let us return to the present case (|I|) of the t-J 
model. The condition that the deconfinement state oc- 
curs may be estimated as (3g'^ < 1. Due to the T- 
dependence of (7|;(T) as shown in ( 2.17 ), the deconfine- 
ment state, hence the CSS, takes place below certain crit- 
ical temperature Tcss; 



C -T^ 
T 



T 



Tc: 



ss 



< 1. 



(2.23) 



For T < Tcss I the MFT works as a first approximation 
because AAxi is small there. We stress that the reason 
why our resu lt T < Tcss for the CSS state is opposite to 
that of ( 2.16| ), T > Tcdpt as in QCD, is du e to t he non- 
trivial T-dependence of coupling constants (2.17), which 
are generated by integrating out holons and spions. In 
Ref.Ej, Nagaosa studied certain dissipative gauge theory 
of fermions and obtained the result that a deconfinement 
state appears above certain temperature or the strength 
of dissipation is strong enough. However, we note that 
he does not discuss the possibility that the parameters of 
his model may be effective and T-dependent. 

At first, the above conclusion may seem rather strange 



6 



since the original gauge coupling is infinite (the coeffi- 
cient of Maxwell term is zero). Our study shows the 
possibility that a coupling to matter-fields strongly in- 
fluences the phase structure. In fact, this fact is now 
well-established in the elementary particle physics, espe- 
cially in the lattice gauge theory. A good example is the 
non-Abelian gauge model without matter fields, which 
is always in the confinement phase in 3D at T = re- 
gardless of the strength ofjlhe gauge couphng. Studies 
of the lattice gauge theoryQ, however, show that, when 
light iVy-flavor quarks are coupled, the phase structure 
drastically changes depending on the value of Nf. Even 
at infinite-gauge-coupling limit, this quark-gluon system 
is in a deconfinement phase for Nf > 7 at T = 0. This 
example is closely related with the present U(l^ gauge- 
theory model for the t-J model. First, PolyakovEZi showed 
that the pure compact U{1) lattice gauge system on a 3D 
lattice without matter fields is always in a confinement 
phase. (This system corresponds to the 2D Hamiltonian 
system at T = 0. In certain approximation, it can be also 
regarded as an effective system of a quantum system in 
3D space at finite T. For more details, see the paragraph 
above (|4.5|).) Next, inclusion of nonrelativistic fermions 
(the spinous in the SB representation) can be regarded 
as inclusion of many-flavor light Dirac fermions because 
of the Fermi surface (line) instead of the "Fermi point" 
of the Dirac fermion. Thus, what wc showed for the t-J 
model in Papers I is nothing but the similar deconfining 
phonomenon induced by ample light fermions. 

Let us make a couple of comments on the explicit curve 
of Tcss('^) in Fig-0- As mentioned in Sect. 2. A, it is much 
lower than the MF value Ty,. This implies that the ef- 
fect of gauge fluctuations are significant. In contrast 
with Fig|lL the observed rcss('5) seems decreasing as S is 
increasedrl In Ref.E3, we observe that the two-body in- 
teractions introduced there tend to lower Tcss for larger 
S. We also note that the value of rcss('5) in Fig.|l] is 
not reliable at very small S , because the long-range AF 
order in 3D at 5 (< 0.04^ upsets the validity of the hop- 
ping expansion we used.B (This point is related with the 
argument given in Sect. 2D below.) 



C. Spin Gap State 

At low temperatures below TsG, the MFT of Sect. 2 
predicts that the NN spin-pairing amplitude A develops 
in addition to x- This phase is called the spin-gap phase, 
because the spinon excitations in 2D acquire the following 
excitation energy Eik), i.e., an energy gap; 




to generate this gap is similar to that of the famous en- 
ergy gap in the BCS model of conventional superconduc- 
tivity. Since the spinon part of the MF Hamiltonian (2.6) 
has the structure, Hf - eflf„ + Kf\fl + fU])^ diag- 
onalization of Hf by a Bogoliubov transformation gives 

rise to Hf ^ Y.k,a E{k)al^aka, E{k) ^ [e^ + X^]^/\ 

In this phase, there app ear t wo kinds of gauge fields 
Axi and B^i of (2.10) and (2.15). Their gauge dynamics 
can be studied by the action Sjjv that is obtained by 
integrating over bx and fxa by the hopping expansion 
w.r.t. Uxi, Vxi- Suv is calculated in Papers I. It takes 
the form; 

Suv — SuVE + SuvM, 

f'^ r 1 

SuVE= drY^ -^2—^f\9rUl,{T)^rUx^{T) 

Jo x% y^EUK^ ) 



5, 



UVM 



•^0 x,t<] 



„2 (rp\UxiUx+i,jUxJ^jiUxj 

i— ^4iV"^^- -Vx+,iV^- 

/rTi\ * XI * x+t 1 X+J,t * XI 



9m V 



-2 77f^X+^J Ux+,^Vx, Ux^ + B.C. 



(2.25) 



Nature of gauge excitations of this system may be 
drawed by calculating the quadratic parts of Sjjvm in 
terms of Axi, Bxi- The g^^UUWW term gives rise 
to the Maxwell term of Axi- If this term alone is 
kept, it describes a massless gauge field Axi- The term 
g'^.^yVV^VV^ by itself also describes another massless 
excitations Bxi- However, the mixed term g^uyV^UVU 
generates the term m\Axi + m%Bxi, which acts as mass 
terms of Axi and Bxi with m\ g cx 'X^X^ ■ Thus the ex- 
citations Axi and Bxi become massive simultaneously at 
T < TsG- This may be called dual Anderson-Higgs mech- 
anism for two kinds of gauge fields; one gauge field acts as 
a Higgs field for the other gauge field and vise versa. So 
we call this phase a Higgs phase, in which the potential 
energy V{xi,X2) between two external charges falls off 
exponentially due to exchange of massive gauge bosons. 
In Table 1, we comp are V (xi,X2) for each phase. The 
lattice gauge model (2.25) is complicated, and no sys- 
tematic study has been done. Instead, in Sect. 4, we shall 
study the nature of phase transition at Tsg in detail us- 
ing the lattice Abelian Higgs model that is derived from 
the effective field theory of Sect. 3. 



D. Lower Dimensional Cases 



We note that there are several gapless points in k space, 
at which cos/ci — cosfc2 — 2jlp / iiJx)- The mechanism 



In Sect.2B, we explained that the CSS can be under- 
stood as a deconfinement phenonenon of the dynamical 
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gauge field A^i- We do not rely explicitly upon the fact 
that the t-J model of our interest has a weak but fi- 
nite three-dimensionality. Then one may ask whether 
the same argument can be applied for systems in lower 
spatial dimensions. We first consider the ID case. Then 
we comment on the 2D case. 

For example, the supersymmetric t-J model (t — J) 
in one-dimension is known to show the CSS.lj However, 
a straightforward application of the arguments in Pa- 
pers I, i.e., the hopping expansion and the Polyakov- 
Susskind theory of CDPTLJ to the ID t-J model leads 
the result that only the confinement phase exists. Tcch- 
nicaly speaking, the Polyakov-Susskind method maps 
the system to an asymmetric XY spin modeO in ID, 
and this spin model is easily shown to have only a dis- 
ordered phase, which corresponds to the confinement 
phaseJj This result seems rather natural because in ID 
the Coulomb force itself generates a linear-rising confin- 
ing potential. There is no possiblity for a deconfincmcnt 
phase to appear in ID if the effective action of the dy- 
namical gauge field is local. However, if the hopping 
expansion is not a good approximation or is not conver- 
gent, a deconfinement phase may exist. This is in sharp 
contrast to the system in spatial 3D in which the exis- 
tence of a deconfinement phase can be predicted in the 
local gauge action and the hopping expansion. 

In order to gain insight into the above problem in ID 
and to show that a deconfinement phase really exists, 
we consider gauge dynamics at low T and move to the 
momentum space instead of the real space assuming the 
fluctuation of the dynamical gauge field is small. This is 
legitimate since at low T the fluctuation of the gauge fleld 
becomes small as in the higher-dimensional cases. Even 
in this case, the gauge interaction gives rise to linear- 
rising conflning Coulomb potential if the effects of mat- 
ters can be ignored. 

We shall focus on the spinon contribution for simplic- 
ity. For Uxi = 1, there exist two Fermi points at ±kp 
whose value depends on the hole concentration. Then 
the left and right moving modes can be defined in the 
usual way. The continuum field theory is easily obtained 
and that is nothing but the 2-flavor massless Schwinger 
model. This model can be exactly solved. Especially 
by integrating out the spinon or the Dirac fermions, the 
effective action of the gauge field is obtained as 



(fx(fx' Af,{x)U^^{x - x')A^{x') (2.26) 



where ^0(2^) = Q^x and Uxi = exp{iAi{x)). The most 
important term in n^i,(fc) is the "mass term" of the gauge 
field. 



H, 



ku 



(2.27) 



The mass term comes from a "massless mode" in the 
particle-anti-particle channel, which is a direct conse- 
quence of the kinematics in ID and behaves like a Higgs 



boson. It is obvious that H^^ (fc) is finite as ^ and so 
the higher-order terms of the hopping expansion in the 
real space are to give the main contribution to it. 

The above consideration in momentum space assuming 
small and smooth fiuctuations of the gauge field gives a 
consistent solution of the deconfinement phase which is 
realized as a Higgs phase. The hopping expansion in the 
real space cannot uncover this possiblity. 

Next, let us consider the 2D case. The method, of Pa- 
per I gives rise to aruasymmetric XY spin modeEldn 2D. 
The MC simulatioicJ and the analytical methocB show 
that this model has an order-disorder phase transition, 
which corresponds to a CDPT. These analyses show that 
the CDPT occurs in higher T than the spin-gap gener- 
ation temperature Tgc as in 3D (See Fig.0). Thus the 
effective model is nothing but the 2D (symmetric) XY 
model, and the CDPT is the Kosterlitz-Thouless transi- 
tion of infinite-order. Although the validity of hopping 
expansion has not yet been clarifled, very existence of 
a CDPT may survive in the correct treatment. There 
is a special point in 2D. The potential energy between 
two slave particles with the opposite charges in Coulomb 
phase is logarithmic V(r) ~ ln(r) due to the 2D Coulomb 
force. Thus the slave particles are "confined" even in 
the Coulomb phase in the sense that limr^oo y{f) = 00. 
This requires a careful study of the nature of quasipar- 
ticles in the Coulomb phase of the pure 2D system. As 
mentioned, inclusion of a finite three dimensionality re- 
leases us from this problem. 



III. EFFECTIVE FIELD THEORY OF THE 
SLAVE-BOSON T-J MODEL 

A. Hamiltonian of the Continuum Field Theory and 
Integration over Holons and Spinons 

Let us consider the low temperature region T < Tqss^ 
where the fluctuations of gauge field Axi is small, as ex- 
plained in Sect.2B. Thus it is reasonable to derive the 
low-energy effective theory there in the form of local 
gauge fleld theory, and apply it the usual perturbation 
theory to evaluate the eff ect of Axi- Let us start with 
the Hamiltonian Hmf of (2.6) and translate the lattice 
variables to the continuum fields such as fxa- ^4 O'faix) 
by using the lattice constant a. For Uxi we setLj 



Uxi = exp(mv4i(a;)) 

~ l-f iaAi{x) - ^a'^Ai{xy 



(3.1) 



The Hamiltonian He of the continuum fleld theory is then 
calculated as 



H,, = I (fx 
(fx 



^^Jx^ + mFx{\Xs{x)\' + \Xd{x)\' 



1 



2mB 



b{x) 
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(fx 



2m p 



lip 



i2ny 



AsG(fc)/|(fc)/|(-fe) + H.c. 



(3.2) 



where (J-b.f denote the chemical potentials for the con- 
tinmnn theory, fa{k) is the Fourier transform of fa{x), 



(3.3) 



and 



\s{x) = -[Xi{x) + X2{x)l 



^d{x) = -^[Xiix) - X2{x)], 



1 



2mB 



2m f 



2tt 



di - iAi, kp = —{I - S) = 2mFfJ- 



(3.4) 



Di is the covariant derivative w.r.t. the gauge field Ai{x), 
and fci? is the Fe rmi momentum of spinons. The last 
term of ffc in (3^) describes the coupling of the spin-gap 
amplitude to spinon pairs. Originally, it has the form of 



/ 



(27r)4 

AsG(fc,q)^2(l-<5) 



AsG(fc,9)/|(fc + |)/|(-fc + |) + H.c. 



2'-'^' 2 



Xdiq)-25XM- (3.5) 



^30(^1 9) consists of the first d-wave term and the sec- 
ond s-wave term. At half-filling, (5 = 0, the s-wave term 
vanishes. In He of (^|^), we simplify Aso by keeping only 
its d-wave term and ignoring its q-dependence as 



AsG(fc) = 2(l-<5) 



kp 



Xri 



(3.6) 



where = Ad(0) ((A^) ~ (— )'Ad). We also set q = in 

n of (H- 

To obtain the effective action Lcs{Ai, Xi) of Xi{x,T) 
and Ai(x, r), we integrate over b and f^ by the standard 
bilinear integrations in the partition function, 



Z = I [db][df][dX^][dA,]e^p{~ J drL,) 

[dXi] [dAi\ exp (- y drLctf {Ai,Xi)), 

Lc = J d^x{b^drb + fldrfa) + H,. (3.7) 



B. Gauge-Field Propagator 

The propagator of the gauge field, Dij{x,T), is gen- 
erated by fluctuations (vacuum polarization) of spinons 
and holons, i.e.. 



D,j{x,t) = {A,{x,t) A, (0,0)) 

= np + iiB] \ 

IlF^Bij{x,T) = ~{Jf,B i{x,T)JF^B j(0, 0))iPi 

+ d,jS{x)d{T)np,B, (3.8) 



where 



l-S 



riB = 



(3.9) 



are the holon and spinon density and 11^,5 ^ represent 
one-particle-irreducible (IPI) diagrams of spinon and 
holon loops made of the following currents coupled to 
Ai ; 



Jfi = Y~ ^i^f'^^ifa + H.C. 



2m p 
1 

2mB 



{ibd,b + B..c.). 



(3.10) 



In the Coulomb gauge, diAi = 0, the propagator at mo- 
mentum q and Matsubara frequency e; = 27r///3 is writ- 
ten as 

A.(g,.O-/(0/fe-^^^^--)A.(..,r) 



D{q,ei) = Ilp{q,ei) +IlB{q,ei) 



(3.11) 



Since we shall need D later in calculating p, we calcu- 
late D below in the random-phase approximation as 



D 



(3.12) 



When the spin gap is sufhciently small, its effect to Hp 
is evaluated by perturbation theory. We obtain 



TT q mp'' ' ^ 



(3.13) 



where we used the relation, fip nnp/mp and vp = 
kp/mp. The second term cx \ei\/q in the first line repre- 
sents the dissipation due to fermions. This term appears 
also in the studies o£,non- Fermi liquids and ^ = 1/2 
quantum Hall effect .ell np{T) is the superfluid density 
of spinons and is calculated for small Asoik) / (kBT) as 



np 
2^ 



AsG(fe) 



2kBT 



ttA 
2kBT 



(3.14) 
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with ki/k2 = tan(/) a nd |fc | = kp- We note that the 
last constant term of ( |3.13 ) represents a mass term of 



Aiix^r). nH is calculated as 



ng(g,Q) 
9^IBi\^iB\) 

fB{e) 



247rm b 



AT) 



, \ei\:>VB{T)q 



exp(/3e) — 1 ' 



"s(^)=T7i h' ^b(T) = ^ . (3.15) 

Eqs.(p3|) and ( p^ ) are obtained for small (/(^ 7r/a). 
For large g's, they should be replaced by anisotropic ex- 
pressions due to AsgC*!)- These anisotropy can be ig- 
nored as long as the spin gap is sufficiently small. 



C. Ginzburg-Landau Theory with Gauge Field 

The most dominant contributions to Z from the in- 
tegrations over Ai and Ai come from their static (r - 
independent) modes. It is seen for Ai by ( 3.13 , 3.15 ). 
So we keep only the static (r-independent) modes in the 
effective Lagrangian density L^s- The partition function 
is then written as 



Z = I Y{Y{d\{x)dA,{x) exp(-/3 ( cfxL^s)- (3.16) 

Loff is given up to the fourth-order in fields and deriva- 
tives by 

= c ^(252|p,.A,|2 + (1 - 5f\V,\d\- 



H.c. 



\2iTm ^ 4 

|2 , „ |\ |2 



+ 26<5^(l-<5)^(4|A.nA,p + A^A^ + A^A^), 

(3.17) 

where the bars represents complex conjugate quantities. 
The GL coefficients, etc. are given by 



1 

m 



1 

nip 



fB{\^iB\) 



2 „ (2e' \ 
mpX 'm-pO In pujx], 



2mB 

2 „o. /2e'^ 

V TT 

mp 2 / 2e'' \ 
ad = mpx (1-d) In pcux] 

TT \ TT / 

b 



7C(3)^2 7C(3) 2,2 
—fi mp, c = oo_.^„^ 



V,, = 



87r3 ' ^ ' " 2,2TT^mp'^ '"^^ 
d^-2iA,, F,.j = d^Aj - djA,, 



(3.18) 



and 7 is the Euler number. uj\ is the cutoff of the spinon 
energy = k^/{2mp) — fip] in the one-loop integrals 
representing spinon pairings, and is estimated as ujx ~ 
Oifip). 

The system favors the d-wave state at small (5's, and 
the s-wave state at large (5's. Actually, the transition 
temperature Tx^^_,^ at which Xd{s) starts to condense is 



estimated from the potential energy V{Xs,\d) of ( |3.17| ) 
by setting ax^f^-^ ~ 0, and one sees that Tx^ > Tx^ [Tx^ < 
TxJ for small (large) (5's. Let us consider the region of 
small (5's and focus on the (i-wave component by setting 
As = and parameterizing as 

Ai(a;) = \exp[i(j){x)], X2{x) = — A exp[i(?!)(a;)]. (3.19) 



Here we introduced the MF A, the (i-wave spin-gap am- 
plitude, for the radial parts of \i{x), ignoring their fluc- 
tuations. Then we have the effective Lagrangian density, 

LcS ^ Lx + La, 



Lx = adX' + (1 - '^)'^'' 



La 



2A, 



127rm ^ 4 



(3.20) 



The MF result is obtained by minimizing the above Lx 
w.r.t. A as 



2e'>' 

^bTsg = exp 

TT 



TTX 



2 \2 



c{i-syx 



k'p 



12TTnip 



(1-5)2 

T \ 



(3.21) 



where Tsg(<5) is the critical temperature below which A 
develops. (See Tsg of Fig.l.) The second result holds 
for T near Tsg- We assume a small but finite three- 
dimensionality to stabilize these MFT results. This point 
and other details of the transition at Tsg are studied in 
Sect.4. 



D. Gauge-Field Mass in Variational Method 

Let us consider the effect of La upon Lx by integrating 
over Ai . We have treated Ai as a noncompact gauge field 
in spite of being originally compact. This procedure is 
appropriate in respect of the kinetic term of Ai , because 
we consider the region T < Tqss- In short, the compact- 
ness of the gauge-field kinetic term should be respected 
in studying the CDPT at Tcss, while the compactness 
of the coupling to the Higgs field is crucial to study the 
transition into Higgs phase at Tsg- So we respect the 
compactness of the mass term by replacing {Ai — di(j>/2)'^ 
in La to the periodic one. Then the new Lagrangian Lb 
which replaces La is expressed as 
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1^1 

127rm ^ A ' 

1-3 



2aB, 



where we introduced the Proca (massive vector) field -B^, 



Bi = A^- -d^cj), 
Fij = diBj — djBi 



(3.23) 



Noe that if the cosine term in Lb is expanded up to the 
second order, La is recovered. Let us take the unitary 
gauge (j) = 0. Then the path integrals reduce as 



lldq> Y[dA, = Y[dq^Y[dB,~^Y[dB, 



(3.24) 



Let us estimate the gauge-field mass tua by the vari- 
ational method.eS We choose the variational Lagrangian 
L'g for Lb as 



(3.25) 



where niA is introduced here as a variational parame- 
ter. There holds the Jensen-Peierl inequality for the free 
energy density F^^'^'^ of Lb, 

^oxact < j^var ^ ^ }_ J ^2^^^^ _ ^,^y ^ 

eM-f3F'B^""'V) = /"j|dB,exp(-/3 f d^xLB), 

qM-PF'bV)^ flldB.expi^-p f d^xL'B), 

{Oy= fY[dB,Oexp(^-P f d^xL'B+pF'eVy 

(3.26) 

where V — J d^x is the system volume. By taking the 
definition of the propagator at the same point and the 
trace of a functional operator O as 



{B,{x)B,{x)) = \iinjB,{x)Bj{y)), 



TrO 



2: Jun (a;|0|y), 



(3.27) 



the variational free energy density Fg^'' is calculated by 
integrating over Bi as follows; 



kBT 
8n 



fcRe = 



1 



= X 



1, 



(3.28) 



where qc is the momentum cutoff of th e Bi modes. Its 
explicit value is estimated later i n (^.21 ) . O is a function 

of ( 3.28| ) we have omitted 
since we are interested 



6(5, T) of S and T. In F^""' 



(3.22) higher-order terms of 0{m\/q'^ 



in small m\/q^. This assumption is justified a posteriori 
by the result itself. 

By minimizing F™{mA) w.r.t. mA, we obtain 



967rm 



d^l 



e 



2d 



2Q-T' 



(3.29) 



We note that the fluctuations of A,- do not affect the MF 



result ( 3.21 ) as long as d > 1 since the corrections then 
becomes higher-order than A"*. Thus the gauge-field mass 
ruA at fixed 5 starts to develop continuously at Tc^g{5) 



as 



m^(5,T)cx(TsG(<5)-r) 



d(S,TsG{S)) 



(3.30) 



That is the exponent d is neither 1/2 nor a constant, and 
drastically changes especially when TsG ~ Ta, where Ta 
is a root of the equation Ta = 6((5, Ta), at which d{5, T) 
diverges. This is in strong contrast with the noncompact 
case, in which the effect of Ai upsets the MF result of GL 
Lageangian Lx as stated. Actually, in the calculation of 
Ref.tj, La in 3D generates a notorious A'^ term. 



IV. PHASE TRANSITION INTO THE SPIN GAP 
STATE: LATTICE ABELIAN HIGGS MODEL 

In this section, we clarify the nature of the phase tran- 
sition at TsG described by the effective field theory (3.2C) 
of Sect.3C. The reader who is interested in the practical 
applications of the effective field theory may skip this 
section, and directly go to Sect. 5. 

To study the transition at TsG, it is useful to intro- 
duce a lattic e mo del that is suggested by the effective 
field theory ( |3.20 ), because ample methods of analysis 
are available for lattice models. The most natural lattice 
model for this purpose is a noncompact version of the so 
called lattice Abelian Higgs model, which is defined by 



/7r i*OQ 
dcf^xW j C^^a;iexp(-S'AH), 

-SaH = -/9^COs[Vi0a; - Axi] + ^ E i^ij^i^x. 



2 

t:3 ] 1 



(4.1) 



where (px € [— tt, tt] is th e pha se of a Higgs field on 



corresponding to 4i{x) of (3.1£) and A^j G [—00,00] is a 



noncompact gauge field corresponding to ^^(a;). The two 
parameters p and k are determined by the parameters 
of the effective model of Sect. 3. (This p should not be 
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confused with the dc resistivity although we use the same 
latter.) ^ah has the local U(l) gauge symmetry under 



A, 



(4.2) 



Its naive continuum limit reduces to ( 3.20| ). 

A similar .but hotter studied model is the compact ver- 
sion of (4.1)c3'La defined by 



/TT PTT 
d^xJl dA„exp(-5iH), 

X x.i ''-^ 



X,l<J 



(4.3) 



It is obtained from (4.1) by (i) letting A 



(ii) replacing the quadratic Maxwell term {eijViAxj) 
by the periodic (compact) one, cos{eijViAxj) 



e [-7r,7rJ, 

^iAxj) 

ViAxj), and (iii) 

introducing q. q is the charge of the Higgs field, and the 
phase structure may differ for each q. On the contrary, 
in the noncompact case of (4.1), we have set q = 1 since 
q can be eliminated by rescaling Ai —> q~^ Ai, k ^ q^n 
and so irrelevant there. 

The main difference between these two models on 3D 
lattice is that monopoles (which are nothing but instan- 
tons in the 3D lattice systerrO) are allowed in ([l.3|), while 
not in ( [4.l[ ) . A monopole is a particular configuration of 

Aa , 



VjAk 



which is a solution of the field equation and has the non- 
vanishing divergence of "magnetic field" Hi. Condensa- 
tion of these monopoles let Ai random, AA^-j^ oo, and 
drive the system into the confinement phase.EZI Since we 
are interested in T around Tsg (< Tcss) which is in 
the deconfinement phase, we consider the noncompact 
model ( [l.lD . On the other hand, we respect the periodic- 
ity (compactness) of the first term in Sah of (^^) , which 
plays the important role in the transition at Tsg- This 
point was already pointed out in Sect. 3D. At the end of 



this section, we shall mention the compact model (4.3). 

Let us first consider the isotropic 3D model (4T) with 
i = 1,2,3, because the phase structure of the quasi-2D 
t-J model in which we are interested is governed by the 
three dimensionality of its effective lattice model. (The 
effect of anisot ropy is discussed later in the paragraph 
containing eq.( 4.14 ) below.) Here we comment on the 
dimensionality of an effective lattice model. It is well 
known that the path integral quantization maps a quan- 
tum system in Z?-dimensional space to a classical system 
in D -|- 1 dimensional space, where the extra one dimen- 
sion represents the imaginary time, t e (0,/3). In many 
cases, this D + 1 classical system is approximated by an 
effective model in Z?-dimensional s pace. A good example 
is nothing but our model, L^g of (3.16), where only the 



T-independent modes appear. As long as T ^ 0, the ex- 
tra dimension has a finite extension /3, and so it should 
not affect the long-range properties and the phase struc- 
ture. (The case T = is different of course.) The T 
dependence are incorporated into the coefficients of the 
effective action. The reduction to the D-dimensional ef- 
fective model is reasonable in this sense. Because the 



Abelian Higgs model is originated from L^s of ( 3.16 ), 
both models should be considered in the same spatial 
dimensionality. 

Let us modify the gauge-Higgs term of S'ah to the fol- 
lowing Villain (periodic Gaussian) form; 



exp 



exp - -{Vi4>j; - Ax^ + 27rn.xiy 



n-ri — — oo 



E 



e"^ > exp 



-TO^i -I- imxi{Vi(i)x - A^i) 



(4.5) 



This replacement keeps the periodicity (compactness) 
in (px- It is then straightforward to perform a duality 
transformationo to reach the following dual representa- 
tion; 

^AH=nE / "^^^ 

xi XI 

(4-4) = Zo n E' ^'^P ( - 4^'pim J-^D,,{x - y- m^)Jyj 



xi J^i xi yj 

(4.6) 

where m? — pj k. The dual variables Jxi (= 0, ±1, ±2, ...) 
and Cxi (g [~oo,oo]) are defined on the link (x,x -\- i) 
of the dual lattice (which is obtained from the original 
lattice by shifting each site x to x + (1 + 2 + 3)/2). From 
the saddle point equations, one obtains the relation. 



(Jxi) oc {eijkVjiVk't'x - Axk)) 



(4.7) 



Thus Jxi represents the gauge- invariant "vortex current" , 
because J (x V x Wcj) for A — 0. (Recall that the vorticity 
in the continuum is defined by f d(f) along a closed line. 
On the lattice, a vortex with a phase change 2-k along a 
plaquette in the 12 plane is expressed by Jx^ = il, and 
so on.) From the gauge inyar.iance, Jxi can be chosen so 
as to conserve at each site, Ho 



^^^Jx 



0. 



(4.8) 



The prime on ^ ^ . of ( |4.6D implies this condition. The 
vector field Cxi mediates the force between two vor- 
tex current elements, Jxi, and generates the potential 



Dij[x — y; m^) as shown in the last line of (4.6) 
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Zq is the partition function of the vector field Cxi with 
mass TO, the Green's function of which is given as 



"^1 ™2 



J 

(-V^ + inn?)D{x - y; m^) = Sxy 



(4.9) 



The system (4.6) is character ized by the density of vor- 
tices. If there are no votices, (4/7) shows that every (px 
is stabilized to a fixed value, say (jjx — 0, with A(j)x — 0, 
and gives rise to a coherent condensation of the Higgs 
field with a long-range order, {Xx) ^ 0. Appearance of 
vortices let (j)x fluctuate, and the prolifiration of vortices, 
i.e., a vortex condensation, leads to destruction of the 
long-range order of the Higgs field. Thus one expects a 
phase transition between these two states; (i) the Higgs 
phase where the Higgs field condenses and vortices do not 
condense, and (ii) the Coulomb phase where vortices con- 
dense and no Higgs condensation, (A^:) = 0. A valance 
between the ene rgy and the entropy determines the crit- 
ical line. From (4.6), the energy of a vortex current with 
the strength \Jxi \ — 1 and the length L is estimated as 



E{L) -47r2pD(0;m2)L, 
whereas its entropy S{L) is estimated as 
S{L) - In/i^, 



(4.10) 



(4.11) 



with some UjUmerical constant /i ~ 4.7 for the simple 
cubic latticccj Thus we estimate the transition temper- 
ature T* (which should be identified with Tsg) as 



F = E{L) -TS{L) c^O, 
T* ~ 47r2pD(0;TO2) In fj,. 

pD{0; TO^) in T* is estimated as 
pD{0;m^) = p [ (fp 



(4.12) 



^ Z^jSin Pi 

p, p/k<1 
K, p/k > 1 ' 



(4.13) 



The duality between the Higgs field X{x) and the vortex 
density v{x) is illustrated in Table 2. 

Detailed studies of the above sjjstem have been per- 
formed by numerical simulationsHj There the radial de- 
grees of freedom of the Higgs field are treated also as dy- 
namical variables. It is found that the phase transiton re- 
ally occurs, and the order of the transition depends on the 
parameters of the model. For a small Higgs-four-point in- 
teraction, it is of first-order, whereas for large foupjDoint 
interaction, it becomes a "continuous" transition.E3 The 
mass of the gauge boson was also calculated in detail 
near the continous phase transition as a function of the 



Higgs mass. In the isotropic 3D case, the gauge boson 
mass develops as (Tsg - T)^^^ as predicted by the MFT 
of (3.21), i.e., TTiA oc A and d — 1/2. 

For the quasi-2D similar phase transiton be- 

tween the Higgs phase and the Coulomb phase is ex- 
pected to occur, i.e., the phase structure of the quasi- 
2D system is governed by its three dimensionality. This 
expectation is quite natural from a general view point 
of universality in renormalization group. The effect of 
the three dimensionality may be described by introduc- 
ing anisotropy parameters in the Hamiltonian such as 



H = - Y^uU+jLfx + .abx+ii.C. 



_ V — I ft - ft ft 

x,i 

ti = t2 = f, fa = at, 
Ji — J2 = J, J3 ~ loiJ. 



(4.14) 



The parameter a controls the quasi-two dimensionality, 
i.e., a = 1 corresponds to the isotropic 3D case and a = 
to the pure 2D case. Another parameter 7 respects 
the difference in ti and Ji. If we use the relation J = 
At^ /U obtained in the derivation of the t-J model from 
the Hubbard model, we have 7 = 0;. 

As an explicit example of the related quasi-2D models, 
we have studied the quasi-2D antiferromagnetic Heisen- 
berg model at finite tJIj By using the quasi-2D 0(3) 
nonlinear tr model, we showed that the critical tempera- 
ture Taf behaves as Taf ~ l/|loga| for small a, which 
vanishes for 2D as expected, and the correlation functions 
exhibit 3D critical behabior in the vicinity of the transi- 
tion point Taf with the interval of AT 0(l/(logQ!)^). 

Yet rSnother example is the the quasi-2D XY spin 
model.c3 There the critical temperature in the 2D case 
remains finite, at which a phase transition of Koterlitz- 
Thouless type takes place. We ipte that the application 
of the Polyakov-Susskind theoryEEl to the effective lattice 
gauge theory of Sect.2B in 2D case maps the system into 
the 2D XY spin model, hence predicts a phase transition 
of Kosterlitz-Thouless type at finite Tsg for 2D. This 
point is different from the Heisenberg model, for which 
the critical temperature vanishes as a ^ 0. 



Returning to the present quasi-2D model (4.1), a phase 
transition should occur at T = Tsg > Oi which may be 
either continuous or discontinuous. In the very narrow 
interval AT around Tsg 7 the system exhibits a genuine 
3D critical properties. This interval, however, may be too 
small to access by the experiments. Outside of this AT, 
the system should exhibit the 2D behaviors, which can be 
studied well by using the results obtained in Sect. 2D. An- 
alytical studies and the MC simulations of the quasi-2D 
case of ( |4.l| ) are welcome to confirm these expectations. 

Let us no w co mment on the compact U(l) Abelian 
Higgs model ( [4.3| ) in 3D. The duality transformation can 
be performed as in the noncompact case.EJ The resul- 
tant partition function is given by the same expression 
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as (^^) but the vortex current Jxi can terminate at lo- 
cations o f in stantons (monopoles) and anti-instantons. 
Thus eq.(O) is modified to 



^ iJxi — QQxi 



(4.15) 



where Qx (= 0,±1,±2, ...) is the instanton density. The 
partition function is obtained by summing over all the 
instanton configurations ui addition to . • Exis- 
tence of sufficient instantons converts the Coulomb phase 
into the confinement phase. 

Recently, this 3D compact U{1) gauge-Higgs model 
was exawiined in the continuum field theory by Nagaosa 
and LeeEa to study the competition of Bose condensa- 
tion versus gauge-field fluctuations. They employed a 
duality transformation similar to the one above on a lat- 
tice. They concluded that, in 3D the vortex condensation 
always takes place for q — 1, and no phase transitions oc- 
cur; the system is always in the confinement phase. Their 
argument does not depend on details of the model and 
the same conclusion is reached for any systems in which 
vortices and instantons exist and a single vortex couples 
to an instanton. However, their conclusion contradicts 
the numerical studiesE3 which exhibit a genuine phase 
transition from the confining phase to the Higgs phase. 
In Appendix, we shall revisit tlais problem and explain 
why the argument given in Ref.E3 is incorrect. There is 
a critical line in the k — p plane, and this line terminates 
at some point. _ 

Finally, we mention the recent work of Balent et al.cj. 
They studied strongly-correlated electron systems in spa- 
tial 2D and argued its phase strucutre by mapping the 
electrons into bosons by using the Chern-Simons gauge 
theory and then by applying the duality transformation 
(See Table 2). They also introduced two composite fields 
for spin and charge degrees of freedom as in the bosoniza- 
tion in ID case. One may expect a similar method may 
be used for the spinous in the SB t-J model. However, 
their way to derive the low-energy theory of the com- 
posite bosons is accompanied with certain assumptions 
and not straightforward. Anyway, we should stress here 
that their approach to the strongly-correlated 2D sys- 
tem is very close to that to the ID system for which 
the bosonization separates the spin and charge degrees 
of freedom quite easily, but loose the ability to describe 
the Higgs phase (i.e., TsG = 0). On the other hand, in 
our approach, we take the three dimensionality of the 
high- Tc cuprates seriously, which is crucial to predict the 
realistic phase structure with Tgc > as we explained 
above. 

In the following section, we shall calculate the conduc- 
tivity in the spin-gap phase. There we assume the quasi- 
two-dimensionality and the genuine phase transition at 
TsG- However, we checked that the kinematical formulae 
appearing in our calculations, which we obtained assum- 
ing TsG > 0, have smooth and weak dependence on a 
near a — ^ From this observation and the estimation 
a ^ 10^^,e3 we set a = in the calculations in Sect. 5 



for simplicity. The errors caused by this simplification 
do not modify the qualitative conclusion there. 



V. RESISTIVITY 

In Sect.5A, we rederive the loffe-Larkin formula for 
conductivity with a general assignement of EM charges 
of holons and spinous. In Sect.5B, the detailed expres- 
sion of the conductivity is obtained. In Sect.SC, the spin 
gap effect on the conductivity is calculated by using the 
gauge-boson mass rriA of Sect. 3D. 



A. Linear Response Theory and loffe-Larkin 
Formula 

Here we shall consider the response of the system to the 
external EM field A^"" {i = 1, 2, 3). We shall work in the 
temporal gauge, where Aq = 0, Aq^ — 0. In principle, 
the effective action of A'l'^ is obtained by integrating out 
all the quantum fields (spinon, holon, dynamical gauge 
filed Ai) and one can calculate the response to A°^ from 
that action. There is an ambiguity in the way how A'^'^ 
couples to spinous and holons. Namely, one may assign 
the EM charge of a holon, Qh and the charge of a spinon, 
Qs arbitrarily as long as the charge of an electon e(< 0) 
is expressed as 



e^Qs- Qhi 



(5.1) 



as suggested by the relation ( |2.2D . r^his problem was 
investigated first by loffe and Larkinell and the result is 
that this ambiguity does not affect the expectation values 
of gauge- invariant physical quantities. For the conductiv- 
ity, this can be seen explicitly in the formula given below. 
Let us recall that the relevant part of the Hamiltonian 
with general charges Qh, Qs is written as 



1 



2mB 
1 



'x}^m^iQhA'i 



iA,)b{x) 



2m F 
{2n) 



^[AsG(fc)/|(fc)/|(-fe) + H.c. 



(5.2) 



Note that the last line acquires no explicit A"^ de- 
pendence since its expression in coordinate space 
~ / dx'^Axif^^fl_^_i remains unchanged under the 
EM gauge transformation, A^i A^j exp(i(6'°'^ -|- 

^x+i))^ fxa ~* fxaGy^vi^QT)- By integrating out the 
spinon and holon fields, we obtain the following effective 
Lagrangian; 

Leff[A„ AH = {A, + Q,AT)n%{A, + Q,Af) 

+ {A, + QhAT)Il%{A, + QhAf) + 0(3), 

(5.3) 
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where 0(3) represents terms higher than quadratic in 
Ai and/or A^^. In order to integrate Ai concretely, 
the above Lagrangian needs to be approximated by a 
quadratic one. However, the 0(3) terms are very impor- 
tant, because they renormahzc the spinon propagator, 
the holon propagator, and the vertices. This effect is 
partially included by replacing H^^^^ with their renor- 

malized quantities, &jl^gy represented by FigjH With 
this replacement, the Lagrangian is expressed as 



Leff[^„ AH = (A, + Q,Ar)n)?(A, 
+ iA + Q,,Ar^)&i{A, 

After integrating over Ai, we obtain 



hAf) 



n = Q^,lls + Q'IIf 



{Qh-Qs? 



e2 ng'+flp' 



(5.4) 



where 11 is nothing but the response function of electrons 

Finally, from the linear-response theory and Leff[^ 
the dc conductivity a{= an — 022) is expressed as 



n 



an = lim lim flniq, — ie), 

So one arrives at the loffc-Larkin formula, 

-1 



a-'=a^' 



a 



F ' 



(5.5) 



(5.6) 



where aB(F) is the conductivity of holons (spinous). 
The last expression of (5.4) is symmetric w.r.t. the 



holon contribution 11^^ and the spinon contribution 11^^, 
and depends only on the difference (5.1). The reason why 
the ambiguity of charge assignment 



Qh 



Qh + c e, Qs Qs + c e, 



(5.7) 



disappears in the final expression of the conductivity is 
that one can make a shift of the integration variable 



A,-c eA'l 



(5.8) 



so as to eliminate the c-dependence from (5.2). By a 
shift of integration variable, the partition function and 
physical quantities are of course unchanged. 



B. Expression of Resistivity 

In the spin-gap state, the spinon conductivity diverges 
CTi? — > 00 due to a superflow generated by the spin-gap 



condensation (A^i) 7^ 0. This is an analog of the well- 
known fact in the BCS theory that the electron resistiv- 
ity vanishes below Tc due to a superflow generated by a 
Cooper-pair condensation. Actually, A^^ has the same 
structure as the BCS model. Thus the total resistivity p 
in the spin-gap state is equal to the resistivity of holons. 



p = a ^ =ag^ 



(5.9) 



Effects of spinous to p is indirect, but certainly exist and 
show up through the dressed propagator D{q,ei) in cal- 
culating lis. 

Now we calculate the response function Hb- By solving 
the Schwinger-Dyson equation approximately according 
to Ref.C3, we obtain 



n_B ii(0, £;) 



E 



d'^q gig J 

(27r)2 mB^ 



Rsiq, £n;£i) 



iei - ieiTBiq, £«; ei) - AEB(g, e„; e;) ' 
Rsiq, e«; q) = Gb((?, en)GB{q, 

„2 



GBiq,£n) = if-n - 



2niB 



(5.10) 



AI]b((7, e„; £;), representing diagrams containing self- 
energy of holons Sb((7, e„), is necessary to keep gauge 
invariance, 

Ai;_B(g, e„;e/) = i?^^(g, e„; e/) x T,Biq,en)G%iq,e„) 

-^B{q,<^n + >^i)G%{q,en + ei) , (5.11) 

However, in the perturbative calculation, this combina- 
tion vanishes in the dc limit by the symmetry under sum- 
mations. We expect this term dose not contribute to the 
dc resistivity, and neglect it hereafter. 

rB(g, e„;e/) represents vertex diagrams, and con- 
tributes to aB- It is given by 



^Biq, e«;e/) 



XruB 



-y 



(2^ 



qx {q' -q) \ q ■ {q' - q) 
\q'-q\ / g2 

X D{\q' - q|,e„' ~ en)RB{q' ,£n';Q), 

(5.12) 

where q x q' = q^q'^ — qyq'x- We fix the length of q in Tb 
to a typical length qB, 



q% = 47rn_B(r). 



(5.13) 



This qB is determined so that the similar integral, ( 5.10| ) 
with an relaxation time in the holon propagator and 
ieiTs -I- ATjb — > 0, being evaluated by this approxima- 
tion, gives rise to the correct result. Furthermore, qB 
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behaves as qs \/2mBkBT at the temperature region, 
ksT ^ nB/iTiB- So this is a natural choice for the typi- 
cal momentu m sca le of holons. We also fix the length of 
q' of D in eg. (5. 12) to (jb to obtain 



rB(gB,e„;ez) 



Sn'^mB P 
X D(qB^{ 
dE - 



sm 



1 - cos^j, £„' - e, 
1 1 



E ien' + iei — E 



(5.14) 



We consider the underdoped region, ns <^ np {5 I), 
and temperatures around (3"^ ~ ub/tub- Assuming that 
D(q,ei) in Tb dominates in the region near the static 
limit ei — 0, we use the upper expressions in eqs.(3.13) 
and ( |3.15 ). In the denominator of D, the dissipation 
term is larger than the q^ term. 



2n\ei\ 



> 



IT q 127rm 



nsiT) 



(5.15) 



as long as ti ^ 0, since their ratio is 
{q%/m)/{y/n/{qB(})) 0[^hB{T)/nF) and smaU. So 
the n'-sum is dominated at e„' = e„. Then we get 



AmriB (3 



sm 



(1 



cos ( 



3mn|(T) 
2mFhB{T) 



2ti 



{sgn(e„ + eO " sgn(e„)} + ©(£?) 



{sgn(e„ + ei) ~ sgn(e„)} (5.16) 



2eiT{T) 



where 



1 



37rm 1 



t{T) 2mB 13 



Smnf{T) 
2'mFfiB{T) 




3rhnf,{T) ^^2 
2mFnB{T) 



(5.17) 



To calculate IIb ij{0, ei) we insert eq.( 5.14 ) into cq.( 5.10| ) 



and do the q-integral and n-sum as in eq.( ^.12 ) to get 
ub iei 



ns»j(0,eO 



TUB iC{T)ei + iT-^{T)sgn{ei)'' 



(5.18) 



where limc,_+o d(T) is finite. Afte r analytic continuation 
> — > ~ie and using eq.( |3.14 ), we finally obtain the 
resistivity as 



rriB 



1 



e^riB t(T) 



cx T 



1 - 



'3/B(|/^B|)m(l-<5) 7r(l-<5)A 



2mi?<5 



2kBT 



(5.19) 



For TsG < T < Tcss, A = Q-4^id this result reproduces 
the T-linear behavior of Rcf.tSlla. For T near and below 
TsG, one expects the following behavior; 



A (X (TsG - T)", 
p^T\l~ciTsG~Tf , 



(5.20) 



where d{6,T) is a "critical exponent" in 2D, which de- 
pends both on T and 6. The downward deviation of p 
from the T-linear behavior is reduced with increasing the 
doping 5. In Fig.|[ we plot p of eq. ([7|) w ith various values 
of d. The MFT value d = 1/2 of ( |3.21| ) is not consistent 
with the experiment which gives rise to smooth deviation 
from the T-linear behavior (no discontinuity in dp/dT at 
T = Tsg). To achieve such a behavior, one needs d > 1. 
This implies that the fluctuation effect of phases of Xxi 
beyond the MFT is important to obtain a realistic curve 

The data of Ref.lij show that one may fit p in a form 
Co -I- CiT for Tsg < T. Xlji^implies spinon contribution 
to p, calculated in Refsl 



{kBT/fiF)^/^ is 

negligibly small compared with cr^^ ^ niBkBT/uB due 
to higher power in T and a small coefficient, aj.^ = for 



T < Tsg as explained, but the discontinuity at Tsg in 
ap^ is not observable due to its smallncss. 

The constant part Co, surviving below Tsg, may be at- 
tributed to scatterings of charged holons with impurities. 
They may be described by -ffimp = ^Vx^l-bx, where Vx 
is a random potential. Actually, standard calculations 
show that it generates T-indepencLeHt contribution to p, 
AcTjj^ oc I/ub at intermediate T'sJla 



C. Spin-Gap Effect on Resistivity 

Let us estimate the momentum cutoff qc of the gauge 
boson. Since the gauge boson is viewed as a composite 
particle made of two spinous and of two holons, it is 
natural to estimate it by using the cutoff kp of spinous 
as 



= 2kf 



(5.21) 



for small S where the spinons dominate over holons. The 
qualitative behavior of p near Tsg is not sensitive to the 
choice of qc- In fact, the exponent c?(Tsg) is independent 
of qc. 

Next, let us consider the renormalization effect of 
the hopping parameter t. We assume that the three- 
dimensional system exhibits Bose condensation at the 
temperature scale of 
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TUB 



(5.22) 



Tb as the observed Tc in the hghtly-doped region. 
Since t - 0.3 eV gives rise to Tb - 3000 K at (5 - 0.15, 
one needs to use an effective hopping parameter t* ^ (LDl 
eV in place of t so as to obtain a reahstic Tc ~ lOOKtil. 

In Fig.^, we present the phase diagram calculated by 
using these parameters, where Tgc is the spin-gap on-set 
temperature, Ta is the temperature at which the mass 
exponent d diverges, and Tb is the bose condensation 
temperature. In Fig.^, we plot the T-dependence of the 
exponent d{T). In Fig.0, we plot p as a function of T 
for several S's. As explained, the curves reproduce the 
experimental data much better than those with c? = 1/2 
of the MF result, showing smooth departures from the T- 
linear curves, i.e., d(TsG) > 1 for the region of interest, 
0.05 <5< 0.15. 



VI. DISCUSSION 

In this paper, we studied the effective gauge theory 
and the phase structure of the slave-boson t-J model. 
We reviewed our analytical studies which indicates the 
CSS at low temperatures as supported by experiments. 
We also explained that the spin-gap state corresponds to 
the Higgs phase in the effective gauge theory. Both ana- 
lytical and numerical studies show that there is a phase 
transition to the Higgs phase in 3D. In the second half of 
this paper, we summerized the calculation of the resistiv- 
ity especially in the spin-gap phase. The results of p{5, T) 
are consistent with the experimental observations. This 
may support that our treatment of gauge-field fluctua- 
tions in a quasi-2D system by the variational treatment 
of compactness is suitable to describe the spin-gap state 
in the t-J model. 

Conccniing to the dimensionality, we cite the work of 
Su et al.E3. They consider the pure 2D SB t-J model, and 
apply i-hfi dually transformation for bosonized spinons as 
in Ref.Ej. They calculate p in the spin-gap state using 
the disorder parameter, the density of vortex condensa- 
tion, which is nonvanishing in the spin-gap state as 
explained in Table 2. (For more details, see ( |A.8| ).) Af- 
ter fitting several parameters to the experimental data, 
they plotted p. We stress that their approach is dual to 
ours, but fails to locate Tsg- This makes their parameter 
fitting rather obscure. 

Finally, we address here on the observability of holons 
and spinons in the CSS state because it often brings some 
misunderstandings. In Sect.2B, we explained that the 
gauge dynamics is in the deconfinement phase or the 
Coulomb phase below Tcss- Then one may ask if the 
holon and spinon appear as quasiparticles, though phys- 
ical quantities are all gauge-invariant. Deconfinement of 
the slave particles does not necessarily imply that these 
particles are "observed" each by each. This point is often 
misunderstood. 



To see this, it is useful to recall QED. In QED in 3-1- 1- 
dimensions, genuine asymptotie states are constructed by 
the gauge-invariant operators^ like 



■(/'^(a;) 2_ \ w{P) exp(ie 



p 



dx^Af^) 



(6.1) 



where V'(a^) is the electron operator and Af^ is the pho- 
ton operator, implies superposition over different 
paths P starting from x to infinity with a suitable weight 
w{P). w{P) characterizes the configuration of electric 
field since the exponential factor acts as a creation op- 
erator of electric flux along the path P. In the decon- 
finement phase like the Coulomb phase, the correlation 
functions of the above operator can be expanded in pow- 
ers of e and at each order of e the infrared singularities 
are concelled out with each other. As the electrons in 
QED, the holons and spinons in the CSS state should ac- 
company clouds of the (soft) gauge bosons and genuine 
quasiparticl es sh ould be described by gauge-invariant op- 
erators like (6.1). However, explicit construction of such 
operators is a future problemEj 

Here we should mention that there is a point that dis- 
tinguishes the quasiparticles of the t-J model from those 
of QED. Since the t-J model is the model of electrons, any 
physical quantities such as conductivity are expressed 
by the electron operators, and therefore gauge- invariant 
(here the gauge invariance means the invariance under a 
local rotation of the phases of slave particles, and nothing 
to do with the EM gauge symmetry). The deconfinement 
of slave particles means that these physical quantities 
that are expressed by slave particles are (approximately) 
calculated by pertubative calcualtion with respect to the 
auxiliary gauge field which couples to the slave particles, 
becuase fluctuations of the gauge field is not large. In 
this sense, the holons, the spinons and the gauge bosons 
can be regarded as (weakly interacting) quasiexcitations 
at low energies that live only inside of the material. This 
point is different from the well known gauge theories 
in particle physics such as QED, in which one can pre- 
pare quasiparticles at low energies as observable incoming 
and outgoing particles at infinitely separated positions in 
space. Thus, we think that one needs further studies to 
propose a clever method to directly measure the charges 
of the slave particles. 

Associated with the above problem of observability of 
the slave partices, the following questiontHi may arise; 
Since the electron operator is constructed as Eq.(2.2) 
one may assign the EM charges of holons and spinons 
freely as long as the electron has the known EM charge 
e(< 0). Because deconfinement implies that the EM 
charges of holons and spinons can be observed separately, 
how does this freedom be fixed ? If this freedom sur- 
vives, it causes nonsensical result like fractional charges 
of holons and spinons. A similar question has been ad- 
dressed in Sect.5A for the conductivity, and the answer 
there was that the conductivity is independent of this 
freedom. 
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Let us consider this question in detail here. This ques- 
tion makes sense if the EM charges, Qh, Qs, which one 
may assign to slave particles, were the absolute charges 
that are measured from the vacuum (the state with no 
electrons). However, the ambiguity in assigning the 
charges of slave particles just reflects our freedom to 
select the reference state from which their charges are 
measured. For example, if we measure the charges of a 
holon, Qh, and a spinon, Qs, from the half- filled state, 
we have Qh = —e,Qs = 0, while if they are measured 
from the vacuum, then Qh = 0, Qs = e. More gener- 
ally, if they are measured from the state of average ab- 
solute charge — c e per site (c is a real number), one has 
Qh — — (— c e) = ce,Qs = e — (— c e) = (1 + c) e. These 



expr essio ns are nothing but what we explained at (5.7) 
and (5.S) in Sect.5A. Thus they are physically equivalent. 
It is also instructive to see that the electron charge per 
site, 



Qc - e ^ Cla, = e ^ flU = e(l - bib.), (6.2) 

a (7 

depends on the state itself. For the half-filled state, 
Yii fi\^) (I*-*) is the state with no slave particles), we have 
Qc = e, while for the vacuum state, ni^!|0)7Qc = 0. 
The EM charge is an additive quantity, and what we actu- 
ally measure as "charges" in the experiments are always 
the differences of the charges. For example, by applying 
the electron variable, Cio- = blfi^, to a state, the EM 
charge changes by Qh — Qs — ce — (1-1- c)e = — e. This 
difference is of course independent of the reference state 
itself and has a physical meaning. 
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APPENDIX A: PHASE STRUCTURE OF THE 
COMPACT ABELIAN HIGGS MODEL 

In Sect. 4, we derived the dual-variable represenation 
(4.6) of the noncompact Abelian Higgs model (4.1). Simi- 
lar analysis is possible in the continuum field theory of the. 
compact Abelian Higgs model (4.3). Nagaosa and LeeO 
studied this model in 3D continuum and concluded that 
the system is always in the confinement phase. In this 
Appendix, we point out that their argument is insufh- 
cient and the phase transition into the deconfining-Higgs 
phase (i.e., the superconducting phase with a massive 
gauge field in their context) certainly exists. 

The action is given by 



S = 



|(V0(x)-A(x) 



^(VxA7 



where we consider the 3D continuum system, and is the 
phase of the Higgs field and A is the vector potential. 
Duality transformation for (A.l) can be performed as 
in the lattice model. We take the U(l) gauge group as 
compact and allow singular configurations of A as 



V • V X A = AtiQ{x), 



(A.2) 



where Q{x) is the instanton (monopole) density as ( 4.15| ) 
in the lattice model. The phase 4'{x) of the Higgs field is 
also compact, so singular configurations with nonvanish- 
ing vorticity are allowed; 



V X V(/i = 27r J. 



(A.3) 



J is the vortex current as (4.7) in the lattice model. If 
the radial degrees of freedom of the Higgs field are intro- 
duced, vortices can appear as regular solitons of the fi^W 
equations, which are called the Nielsen-Olesen vortexa. 
In any way, it is standard to introduce a complex scalar 
"vortex field" ipv, whose world lines are identified with J. 
The field il)v{x) is viewed as the creation(annihilation) 
operator of a vortex. Then il)y{x)^liv{x) measures the 
vortex density. V'v interacts with each other via a short- 
range interaction as in (4.6). 

One can sum up instantons by the dilute gas approx- 
imation. To do this, one introduces a phase field ^p{x). 
In fact, the relation (A.2) is respected by the Lagrange 
multiplier field <p(x) as 



d(y3(a;)exp(^i(p(x)(^V • V x A{x) ~Q{x))], (A.4) 



and the summation over the instanton Q{x) = 0, ±1 gives 
rise to a factor, 



E 

Q(a:)=0,±l 



l + z 
exp 



(e-*'^("Vy(a;) + e''^(")i^t,(x) 



(A.5) 



where z is the fugacity of instantons, and the first line 
describes that each instanton Q{x) is connected with 
the vortex current generated by ipvix). (V'y (a;)'^^^-' = 
(-0y(a;)'l')~'5('") for Q{x) < 0.) As the factor exp{iLpQ) in- 
dicates, the angle ip{x) is conjugate to the integer Q{x). 
Thus, there holds the uncertainity principle, 

A^(a;)AQ(a;) > 1. (A.6) 

This implies that the following two phases are possible; 

(i) Confinement : Aip{x) = 0; AQ{x) — oo, 

(ii) Deconfinement : A(p{x) = oo; AQ{x) = 0. (A. 7) 



(A 1) (i) is the confinement phase where the instanton con- 
denses, hence Ai fluctuates violently AAi = oo as seen 
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by (A. 2). (ii) is the deconfinement phase where no in- 
stantons appear Q{x) = 0, AQ{x) = 0, so AAi = 0. 

Next, we note that the Higgs field exp(«0) and the 
vortex field ipv are "conjugate" each other in the sense 
that the following two cases are possible; 

(a) Disorder — Vortex : {tpv) ^ 0; {exp{i(j))) — 0, 

(b) Order - Higgs : (V-y) = 0; (exp(i(/))) ^ 0. (A.8) 

Actually, when vortices proliferate with finite density 
{iPy{x)iPv{x)), {ipv{x)) ^ 0, while — oo and so 
(exp(i0)) = 0. This is the case (a), which one may call 
the disorder phase, Acj) = oo, or the vortex phase. On the 
other hand, when vortices are rare, (tpvix)) = 0, while 
(exp(i(/))) 7^ due to At/) = 0. This is the case (b), which 
on may call the ordered phase, A0 = 0, or the Higgs 
phase. This Higgs- vortex duality is summerized in Table 
2, where the vortx density v{x) — {ipY{x)ipvix)) ■ 

The continuum action can be ejcpressed in terms of ip 
and the dual variables as foUowstJ 



d-^x 



1 



(V X 



1 

2k 



1 

2k 



(A.9) 



By integrating over ^jy, its world lines (the vortex cur- 
rent) are generated, and with each line segment the cur- 
rent c is coupled. Integration over the force-mediating 
field c produces interactions among the vortex cur- 
rents. The mass of the vortex field is given as = 
Att^ pD{Q;m^) — In/i as in the lattice model. There is 
also a short-range repulsive interaction (g > 0) between 
vortices. 

To study the phase srucutre, Nagaosa and LeeE3 intro- 
duced a new field ip'y = V'ye"''^ by a field redefinition 
or a "gauge transformation" . The last line of the action 



( [A.9| ) then generates a linear term of ipy. Due to this 
term, they simply concluded that ^p'y always condenses 



(A.IO) 



regardless of the value of A/^. From this, they argued 
that the system is always in the vortex condensed phase 
of (a) of ( [A.8| ) , and there occurs no phase transition into 
the Higgs phase (b) of (^J). 

However, their argument above using the field redef- 
inition or gauge transformation is not correct. A good 
counter example is the gauge-Higgs model itself. In the 
original variables, the gauge-Higgs coupling is given by 



gi0(a;4-i)g«'4i(3:)g- 



-fH.c. 



(A.ll) 



We can fix the gauge to the unitary gauge, i.e., 
exp(z0(a;)) — 1 and then the linear term of the expo- 
nentiated gauge field exp(zAi(x)) appears in the action, 
i.e., S ~ UUUU + U. But this linear term does not 



necessarily lead to {exp{iAi{x))) ^ which is the condi- 
tion that the system is in a nonconfining phase, like the 
Coulomb phase or the Higgs phase. For example, a MFT 
of this model with a variational action of single-link form 
S ^ U offers us two possibilities,^ (i) (U) ^ for the 
Coulomb or the Higgs phase or (ii) (U) = for the con- 
finement phase, and predicts a phase transition between 

(i) and (u). 

To confirm this point, let us analyse the system ( |A.9| ) 
in a mjMe straightforward manner referring to Polyakov's 
result E3 He considered the dynamics of the pure compact 
U(l) gauge model in 3D[{2 + 1)D], and showed that the 
system is always in the confinement phase. After integra- 
tong over the gauge field (instantons) , the action <S'c/(i) 
is given by 



r 1 , - 



2k 



z cos ip 



(A.12) 



where the field tp here again conjugates to instanton 
density and mediates the force among instantons. One 
can approximate the relevant potential term cos ip as 
cos ~ to find that the vacuum is given hy p> = 
(mod 27t). This obviously means A(p <C 1, and so a 



condensation of instantons; the state (i) of (A. 7) above. 
Thus the system iSJiU the confinement phase. This is the 
result of Polyalov.E3 



Now let us return to the present model (A.9) and con- 
sider two cases, JVP < and > 0, separately. We 
assume z is sufficiently small, so that the renormaliza- 
tion effect of Af^ is negligible, although the last z term 
gives rise to a negative renormalization to Af^. 

For AP < 0, the vortex field condenses {ipy) = Vo 7^ 
because the effective action of ipv has a negative curva- 
ture at the origin as in the standard Ginzburg-Landau 
theory with the global U(l) symmetry. Then the field ip 
acquires a potential term like zipo cos ip where we assumed 
that tpo is real without loss of generality. The instanton 
action reduces to that of the pure U(l) gauge model 



( A.12 ), and the Polyakov's result applies; condensatioa 
of instantons occurs and the confining phase is realized.!^ 
This phase corresponds to the confinement-vortex phase 



(ia) of (|Aj) and (|A. 

For Af^> we can safely integrate out the vortex field 
because there are only short-range interactions between 
vortices. Then the resultant action describes the dynam- 
ics of instantons or of their conjugate variable ip{x). The 
action of ip{x) is given as 



2k 



(A.13) 



The second term in (A.13) originates from the fact that 



an instanton and an an ti-instanton is connected by a vor- 
tex flux tube via ( 4.15 ). We recognize that, for large Af , 
the second term of (A.13) is approximated by a derivative 
term, hence 
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(A.14) 



There is no potential term of in (A.14) and Lp does 
not have a definite value, but fluctuates randomly (for 
sufficiently large k and small z^). This is the state (ii) 
above, leading to a deconfinement phase. Here, we have 
{ipv) = because the large fluctuation of p makes the 
last linear term in (pvix) irrelevant; (exp(i(^)) = 0. From 
( [A.8| ), this means that the original Higgs field has a coher- 
ent phase. This phase corresponds to the deconfinement- 
Higgs phase (lib) of ( |A7| ) and ( |A^ ). 

The above arguments indicate that there is a phase 
transition from the confining to the Higgs phases. This 
is actually found now in the numerical studies.Ej On the 
other hand, for the compact gauge-Higgs model with 
fixed radius of the Higgs field, Fradkin and Shenker 
showecEj that the confinement and Higgs phases are con- 
nected if the Higgs charge is fundamental q — ±.\. The 
above two observations can be reconciled if the critical 
line terminates at some point in the (k — p) plane. If the 
radius of the Higgs field is also a dynamical variable, the 
critical line does not terminate as the parameter space is 
now three dimensional instead of two. 
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Phase 
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V{r) 


Confinement 


Tcss < T 


oc r 


Coulomb 


TsG <T < Tcss 


r 


Higgs 


T<TsG 


OC exp(— TO^ r)/r 
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Table 1. Phases of the gauge dynamics of the SB t-J 
model. V(r) is the potential energy between two external 
charges separated by a distance r. 




(b) 



Phase 


T 


X{x) 


v{x) 


Coulomb 


TsG <T < Tcss 





^0 


Higgs 


T<TsG 








Table 2. Comparison of the ordinary represenation 
and the dual representation of the 3D (quasi-2D) non- 
compact Abelian Higgs model for the transition at Tsq 
into the spin-gap state. The order parameter of the ordi- 
nary representation is the spin-gap amplitude A(x), while 
the disorder parameter in the dual representation is the 
vortex density v{x). In the pure 2D system, Tsg = and 
only the Coulomb phase exists, in which v(x) is always 
nonvanishing. 



FIG. 2. Illustration of electric fluxes connecting two ex- 
ternal charges, (a) Conflnement state, (b) Deconflnement 
state. 
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FIG. 3. Graphical representation of the Schwinger-Dyson 
equation for Ilp^^gy 
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FIG. 1. Phase diagram of the SB t-J model in the S — T 

plane. Along and Trvb, X ^i^d D vanish, respectively. 
Tbc is the onset T of Bose condensation, {b^} / 0. Tcss 
is the critical temperature below which the CSS takes place. 
There are five phases: (i) Strange Metal Phase: Deconfine- 
ment-Coulomb; (ii) Spin-Gap Phase: Deconfinement-Higgs 
with D 7^ 0; (iii) Fermi-Liquid Phase: Deconfinement-Higgs 
with (bx) 7^ 0; (iv) Superconducting Phase: Deconfine- 
ment-Higgs with D, (bx) 7^ 0; (v) Electron Phase: Confine- 
ment above Tcss- 
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FIG. 4. Plot 

of the resistivity p divided by po = 2TimFkBTsG / {e^ np) for 
\s{T) ~ Ao(f - T/Tscf (d = 1/2, 1,2). For definiteness we 
chose 5 = 0.05, Ao = IksT* j-K, 2-KnB/mB = ASksTsG/ {l~ 5) 
and J/t = 0.35. 
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FIG. 5. Mean-field phase diagram of the t-J modeL Tbc 
is the Bose condensation temperature. T\ is the spin-gap 
onset temperature. Ta is the root oi T = Q{T) at which 
d{T) diverges. We chose t* = 0.01 eV, J = 0.15 eV, and 
= 7rJx/(2e^). 
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FIG. 6. Exponent d{T) vs T for 5 = 0.04 and 0.1. 
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FIG. 7. Resistivity p{5,T) in h/e'^ for several S's with the 
parameters chosen in Fig. 2. The dotted lines represent the 
case of X{T) — in (1). The exponent d(rsG) decreases as 
16.4, 4.8, 2.8, 2.0, as 5 increases. 
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